ON THE REALITY OF THE EIGENVALUES FOR A CLASS OF 
PT- SYMMETRIC OSCILLATORS 



K. C. SHIN 

Abstract. We study the eigenvalue problem —u"{z) — [{iz) m + P(iz)]u(z) = Xu(z) with 
the boundary conditions that u{z) decays to zero as z tends to infinity along the rays 
argz = — ^ ± ^^2, where P(z) = a\z m ~ x + aiz m ~ 2 + • • • + a m _iz is a real polynomial and 
m > 2. We prove that if for some 1 < j < ^, we have (j — k)a,k > for all 1 < fc < m — 1, 
then the eigenvalues are all positive real. We then sharpen this to a slightly larger class of 
polynomial potentials. 

In particular, this implies that the eigenvalues are all positive real for the potentials 
aiz 3 + (3z 2 + ~fiz when a, /3, 7 € M with a ^ and a 7 > 0, and with the boundary 
conditions that u{z) decays to zero as z tends to infinity along the positive and negative 
real axes. This verifies a conjecture of Bessis and Zinn- Justin. 

Preprint. 

1. Introduction 
1.1. The main results. We are considering the eigenvalue problem 

(1) -u"(z) - [{iz) m + P{iz)\u{z) = Xu{z) 

with the boundary conditions that u(z) decays to zero as z tends to infinity along the rays 
argz = -f ± where m > 2, A G C and P is a real polynomial of the form 

(2) P(z) = a x z m - x + a 2 z m - 2 + ■■■ + a^z, with all a k G R. 

The boundary conditions here are those considered by Bender and Boettcher |IJ. Note that 
the boundary conditions for m = 3 are equivalent to the conditions that u decays to zero as 
z tends to infinity along the positive and negative real axes. If a non-constant function u 
along with a complex number A solves ([!]) with the boundary conditions, then we call u an 
eigenfunction and A an eigenvalue. 

Before we state our main theorem, we first introduce some known facts by Sibuya ]19[ 
about the eigenvalues A of ([!]), facts that hold even when a k G C. 

Proposition 1. The eigenvalues A& of (|7[) have the following properties. 

(I) Eigenvalues are discrete. 
(II) All eigenvalues are simple. 
(Ill) Infinitely many eigenvalues exist. 
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(IV) Eigenvalues have the following asymptotic expression 
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where the error term o(l) could be complex. 

We will give precise references for Proposition [j] after Proposition ^ in Section |2|. In this 
paper, we will prove the following theorem that says that the equation (|1]) with a polynomial 
potential in a certain class has positive real eigenvalues only 

Theorem 2. Let a k s be the coefficients of the real polynomial P(z) = aiz" 1 ^ 1 + a 2 z m ~ 2 + 
• ■ • + a m _iz. If for some 1 < j < y ; we have (j — k)a k > for all k, then the eigenvalues 
of (|7|) are all positive real. 

Corollary 3. In particular, with m = 3 the eigenvalues A of 

-u"{z) + (iz 3 + (3z 2 + 'yiz)u(z) = Xu(z), u(±oo + Oi) = 0, 
are all positive real, provided f3 £ M and 7 > 0. 

Proof. This is a special case of Theorem ^ with m = 3, j = 1 and -P(-z) = /?£ 2 — 72. □ 

We also mention that Delabaere et al. || ||| studied the potential iz 3 + r yiz and showed that a 
pair of non-real eigenvalues develops for large negative 7. And Handy et al. []13 , |14| showed 



that the same potential admits a pair of non-real eigenvalues for small negative values of 
7 w -3.0. 

Remark. By rescaling, the conclusion of Corollary |] holds for the potential aiz 3 + [3z 2 + / yiz 
when «6l- {0}, (3 E R and 07 > 0. 

This paper is organized as follows. In the rest of Introduction, we will briefly mention 
some earlier work. Then in the next section, we state some known facts about the equation 
(H) and examine further properties. In Section ^ we prove Theorem |2], and in Section § we 
extend Theorem ^|. Finally, in the last section we discuss some open problems for further 
research. 

1.2. Motivation and earlier work. Around 1995, Bessis and Zinn- Justin conjectured that 
eigenvalues of 

f_ 

dz 2 

are all positive real. And later Bender and Boettcher |IJ generalized the BZJ conjecture; 
that is, they argued that eigenvalues of 

d 2 
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are all positive real when f3 > 0. Notice this follows for (3 < by Theorem |2| with P(z) = f3z 2 . 
The case (3 > is open, except for m = 3, 4 which are covered by Theorem 0. 
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Recently, Dorey et al. [10], |Tl| have studied the following problem 

d 2 



dz 2 
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with the boundary conditions same as those of ([I]), and M, a, I being all real. They proved 



2/+l|, 
= 0), 



that for M > 1, a < M+l+|2/+l|, eigenvalues are all real, and for M > 1, a < M+l- 
they are all positive. A special case of (^) is the potential iz 3 (when M — |, a = 
which is the (3 = version of the BZJ conjecture, but their results do not cover the (3 ^ 
version. (Suzuki |20) also studied the whole I = version of (0) under different boundary 
conditions.) 

The proof of our main theorem, Theorem El has two parts. The first part follows closely 



the method of Dorey et al. pT|, |r2|, developing functional equations for spectral determinants, 
expressing them in factorized forms and then studying an "associated" eigenvalue problem. 
We also introduce a symmetry lemma that is required by our more complicated potentials. 
The second part builds on earlier work of the author in estimating eigenvalues of the 
"associated" problem by integrating over suitably chosen half-lines in the complex plane. Of 
course both this paper and [JTlj, |12| are indebted to the work of Sibuya . 

Note that our result Corollary ||] proves the full BZJ conjecture; that is, eigenvalues A of (f|) 
are all positive real. Also Theorem |2| contains the polynomial potential case (I — 0, M 6 N) 
of problem @, though only with a < 0, whereas Dorey et al. handle a < M. (Our proof in 
the case a < can be seen to reduce to that of Dorey et al. ) In Theorem [IT] we do manage 
to handle the case < a < M, by using also the harmonic oscillator inequality, which is a 



different approach from that used in [11 , page 5701]. 

In a related direction, Bender and Boettcher 0] found a family of the following quasi- 
exactly solvable quartic potential problems 
d 2 



(7) 



dz 2 



[(iz) 4 + 2a(iz) 3 + (a 2 - 2(3){iz) 2 - 2(af3 - J)(iz)] 



ulz 



Xu(z) 



with the same boundary conditions as those of ([!]), where a, (3 £ R and J G N. Note 
here that the positive integer J denotes the number of the eigenfunctions that can be found 
exactly in closed form. However, for the purpose of studying the reality of the eigenvalues, 
we can allow J 6 R. Our results in Theorem ^| confirm that if for any J £ R, we have either 
a (3 > J and a > 0, or a (3 > J and 2(3 > a 2 , then eigenvalues of (0) are all positive real. 

The above Hamiltonians are not Hermitian in general. However, according to Bender 
and Weniger ||, Hermiticity of traditional Hamiltonians is a useful mathematical constraint 
rather than a physical requirement, in order to guarantee real eigenvalues. All Hamilto- 
nians mentioned above are the so-called PT-symmetric Hamiltonians. A PT-symmetric 



4 



Hamiltonian is a Hamiltonian which is invariant under the product of the parity operation 
V{: z i— > —~z) (an upper bar denotes the complex conjugate) and the time reversal operation 
T(: i i — > —i). These PT-symmetric Hamiltonians have arisen in recent years in a number 
of physics papers, see |7], [13|, [14|, [16], [17|, |21|] and other references mentioned above, which 
support that some PT-symmetric Hamiltonians have real eigenvalues only. In general the 
PT-symmetric Hamiltonians are not Hermitian and hence the reality of eigenvalues is not 
obviously guaranteed. But the important work of Dorey et al. |TTJ , and results in this paper, 
prove rigorously that some PT-symmetric Hamiltonians indeed have real eigenvalues only. 

As a final remark of the introduction, we mention that if H = — + V(z) is VT- 
symmetric, then V(—~z) = V(z) and so Re V(z) is an even function and ImV(z) is an odd 
function. Hence if V(z) is a polynomial, then V(z) = Q(iz) for some real polynomial Q. 
Certainly (|l|) is a PT-symmetric Hamiltonian. 



2. Properties of the solutions 

In this section we will introduce some definitions and known facts related with the equation 
([!]). One of our main tasks is to identify the eigenvalues as being the zeros of a certain entire 
function, in Lemma || But first, we rotate the equation ([!]) as follows because some known 
facts, which are related to our argument throughout, are directly available for this rotated 
equation. 

Let u be a solution of (jl]) and let v(z) = u(—iz). Then v solves 

(8) -v"{z) + [z m + P{z) + AM*) = 0, 

where m > 2 and P is a real polynomial (possibly, P = 0) of the form 

P(z) = a x z m - 1 + a 2 z m ~ 2 + ■■■ + a m _iz. 

Next we will rotate the boundary conditions. We state them in a more general context by 
using the following. 

Definition. The Stokes sectors Sk of the equation (H) are 

f 2/C7T 7T 1 

Sk = s z e C : I arg z 1 < > for k E Z. 

[ m + 2 m + 2J 

See Figure [EJ 

It is known that every non-constant solution of (§) either decays to zero or blows up 
exponentially, in each Stokes sector Sk- Thus the boundary conditions on u in (|l|) become 
that v decays in S_i U Si. 
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Figure 1. The Stokes sectors for m = 3. The dotted rays are argz = ±|, n. 

Before we introduce Sibuya's results, we define a sequence of complex numbers bj in terms 
of the afc and A, as follows. For A G C fixed, we expand 

(1 + a x z~ x + a 2 z~ 2 + ■■■ + a m . lZ x - m + \z~ m ) 1/2 

= 1 + {aiz- 1 + a 2 z~ 2 + ■■■ + a m ^z x - m + Xz~ m ) k 

k=i ^ ' 

°° b 

(9) = 1 + for large \z\. 

i=i 

Note that bi, b 2 , ■ ■ ■ , 6 m -i do not depend on A. We further define r m = — ^ if m is odd, 
and r m = — ^ — bm. + i if m is even. 

Now we are ready to introduce some existence results and asymptotic estimates of Sibuya 
TJ| . The existence of an entire solution with a specified asymptotic representation for fixed 



afc's and A, is presented as well as an asymptotic expression of the value of the solution at 
z = as A tends to infinity. These results are in Theorems 6.1, 7.2 and 19.1 of Sibuya's book 
T9| . The following is a special case of these theorems that is enough for our argument later. 



The coefficient vector 

a := (d, a>2, . . . , a m -i) 

is allowed to be complex, here. 

Proposition 4. The equation (^), with G C, k = 1, 2, . . . , m — 1, admits a solution 
f(z, a, A) with the following properties. 

(i) f(z, a, A) is an entire function of (z, a, A). 



(ii) f(z,a,X) and f'(z,a, A) = T~f(z, a, A) admit the following asymptotic expressions. Let 
e > 0. T/ien 

/(z,a,A) = z^(l + 0(z-^ 2 ))exp[-F(z,a,X)], 

/'(*,a,A) = -^ + ^(l + 0(2- 1 / 2 ))exp[-F(^a,A)], 

as z tends to infinity in the sector \ &rgz\ < — e, uniformly on each compact set of 



(a, X) -values . Here 
F(z, a, A) 
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(iii) Properties (i) an<i (ii) uniquely determine the solution f(z, a, A) o/ (Qj. 

(iv) For each fixed a and 5 > 0, f and f also admit the asymptotic expressions, 

(10) /(0,a,A) = [l + o(l)]A- 1 / 4 exp[irA5+^(l+o(l))" , 

(11) f(0,a,A) = -[l + o(l)]A 1 / 4 exp[irA^ + ^(l + o(l))" 
as A tends to infinity in the sector \ argA| < tt — 5, where 

(12) K = So ( v/1 + tm ~ ^) dt - 

Proof. In Sibuya's book [|I^], see Theorem 6.1 for a proof of (i) and (ii); Theorem 7.2 for a 
proof of (iii); and Theorem 19.1 for a proof of (iv). And note that properties (i), (ii) and 
(iv) are summarized on pages 112-113 of Sibuya's book. □ 

We now give references for the proof of Proposition |l|. We use the number 

2tt 



uo = exp 



m + 2 



(13) 



Proof of Proposition^. See Theorem 29.1 of Sibuya [19] for a proof which says that eigen- 
values are simple, and 

2m 

'-2k + l)vr\ ™+ 2 
2K sin — / 

where K is given by (|12|) . Note that Sibuya studies the equation (||) with the boundary 
conditions that v decays in SoU^, while in this paper we consider the boundary conditions 
of the rotated equation (||]) that v decays in SLi U S±. The factor u m in our formula (13) is 
due to this rotation of the problem. 

The remaining two claims (I) and (III) are easy consequences of the asymptotic expression 
(0)- 



Also one can compute K directly or see the equation (2.22) in [12], which says 



K 



2^ 



r 



m 



r i + - 



So this along with (|T3| ) and the identity r(A)T(l — A) = 7rcsc(7rA) implies (0). Note that 
the asymptotic expression (|3|) of the eigenvalues agrees with that of Bender and Boettcher 
IJ obtained by the WKB calculation for the eigenvalue problem (|5|), after an index shift. 

We mention that the simplicity of the eigenvalues can be proved by using the fact that 
for each Stokes sector, there exist two solutions of (§) with no boundary conditions imposed 
such that one decays to zero and another blows up as z tends to infinity in the sector. □ 

The next thing we want to introduce is the Stokes multiplier. First, we let 

G k (a) := {uj~ k ai , uj- 2k a 2 , ... , (m ~ 1)fc a m _i) for k e Z. 

Let f(z, a, A) be the function in Proposition |]. Note that f(z, a, A) decays to zero expo- 
nentially as z — > oo in So and so f(z, a, A) blows up in S-i U Si. Then one can see that the 
function 

f k (z,a,X) := f(u- k z,G k (a),u- mk X), 
which is obtained by rotating f(z,G k (a),uj~ rnk X), solves (^). It is also clear that f k (z,a,X) 
decays in S k and blows up in Sk-i U S k +i since f(z, G k (a) , uj~ mk A) decays in Sq. Then since 
no non-constant solution decays in two consecutive Stokes sectors, f k and f k +i are linearly 
independent and hence any solution of (|8[) can be expressed as a linear combination of these 
two. Especially, for some coefficients C(a, A) and C(a, A), 

(14) f-i{z, a, A) = C(a, X)f {z, a, A) + C(a, \)fi{z, a, A). 

These C(a, A) and C(a, A) are called the Stokes multipliers of /_i with respect to f and fi. 
We then see that 

W_i,i(a, A) ~ ty_i,o(a, A) 

W 0A (a, A) W 0)1 (a,X) 

where Wj^ = fjf'k ~ fjfk is the Wronskian of fj and fk- Since both fj, f k are solutions 

of the same linear equation ([8p, we know that the Wronskians are constant functions of z. 

Since fk and fk+i are linearly independent, Wk,k+i 7^ for all fceZ. Moreover, we have the 

following which is needed in the proof of our main theorem. 

Lemma 5. The Stokes multiplier C (a, A) is independent of X. Moreover, if a e M m_1 then 
\C(a,X)\ = l. 

Proof. First note that Sibuya's multiplier c(a, A) in ]T9| is Wip/Wip while we use C(a, A) = 
W-i )0 /W ,i- Since f k (z, a, A) = f(u~ k z,G k (a),LU~ mk X), we see that 

fk+i(z, a, A) = /(c-( fc+1 )2,G fc+1 (a),^ fc+1 )A) 
= f k {u- l z,G{a),u>- m X). 

Hence using uo m+2 = 1, we see that 

(15) W k+ i, j+1 (a,X) =uj- 1 W k , J (G(a),uj 2 X), 



which is the equation (26.28) of |T{|. 

So using the equation (26.29) on page 117 of fl9j , one can get 

W-i,o(M) 



C(a, A) 



W ,i(a,X) 

W ,i(G'- 1 (a) ) a;- 2 A) 



" W h2 {G-\a)^\y by(15) ' 

(16) = -^-^(G-Ha))^ by (26>2g) of p 



where 



~(r- l ( \\ — / ^ if 771 is odd, 

^ \ frm/2+i(G'~ 1 (a)) if m is even. 



From fll6|) , it is clear that C(a, A) is independent of A. We want v[G 1 (a)) to be real if 
a G R m -\ so that |C(a, A)| = \u\ = 1. Suppose a G R m_1 . Since fem+^G-^a)) = -6m +1 (a) 
as noted on page 117 of |L9| (or can be directly verified from (||)), it is sufficient to show that 
bm. + i(a) is real when a G IR m_1 . Since a^'s are all real, from we conclude that &™ + i(a) 
must be real. This completes the proof. □ 

Thus from the proof of Lemma [5] we get C(a, A) = for some 0o = 0o(&) G M and hence 
from flli]) we have 

(17) C(a,A)/ (z,a,A) = a, A) - e^f^z, a, A) 

(18) = f(wz, G~ 1 (a),u m X) - e l ^ Q f{uj~ l z, G(a),u~ m X). 

From this, for each a G M m_1 we can relate the zeros of C(a, A) with the eigenvalues of (|I]) 
as follows. 

Lemma 6. For each fixed a = (a±, a-i-, ■ ■ ■ , a m -i) G IR" 1-1 , a complex number X is an eigen- 
value of (flp if and only if X is a zero of the entire function C(a, A). 

Hence, the eigenvalues are discrete because they are zeros of a non-constant entire function. 
Note that the Stokes multiplier C(a, A) is called a spectral determinant or an Evans function, 
because its zeros are all eigenvalues of an eigenvalue problem. 

Proof. Suppose that A is an eigenvalue of (|T]) with the corresponding eigenfunction u. Then 
we let v (z) = u(—iz), and hence v solves (B) and decays in S-x U Si. Since /_i is another 
solution of (H) that decays in SLi, we see that f_i is a multiple of v. Similarly fi is a multiple 
of v. Hence the right-hand side of (|17|) decays in S-i USi. But / blows up in 5_i U Si, and 
so (|17D implies C(a, A) = 0. 

Conversely we suppose that C(a, A) = for some A G C. Then from flT7|) we see that /_i is 
a constant multiple of /i. Thus both are decaying in 5_i U 5i and hence u(z) := f_i(iz, a, A) 
is an eigenfunction of ([!]) with the corresponding eigenvalue A. □ 



Next we examine fll8f) and its differentiated form at z — 0, which are, 

(19) C(a,A)/(0,a,A) = /(0, CT^a), a; m A) - e^/(0, G(a), ^" m A), 

(20) C(a, X)f'(0, a, A) = w/'(0, G _1 (a), cu m A) - e^uT 1 /' (0, G(a), uT m A). 

The right-hand sides of these are given by differences of two functions of A. We will express 
these right-hand sides with single functions, respectively. To this end, we prove that / and 
/' both have some symmetry as follows. 

Lemma 7. Let a = (ai, a 2l . . . , a m -i) G C m_1 . Then we have 



(21) /(0,a,A)=/(0,a,A) and /'(0, a, A) = /' (0, a, A). 

Especially, we have that if a = (a 1; 02, . . . , a m _i) G R m_1 zs rea/ ; i/ien 



(22) /(0,G(a),A) = /(0,G- 1 (a),A) and /'(0, G(a), A) = f'(0, G _1 (a), A). 

Proof. Let g>(z) = /(z, a, A), which is the entire function / in Proposition [|and hence decays 
in 5*0- Then g solves 

-g"{z) + {z m + a x z m ~ x + a 2 z m - 2 + ■■■ + a m _ x z + X)g(z) = 0. 



Next we take the complex conjugate of this and replace z by z. Then we see that g(z) is 
entire and solves the following equation 



(23) -g"(z) + (z m + a{z m - v + ^ m " 2 + • • • + a^T x z + X)g(z) = 0. 



Since the entire functions g(z) and f(z,a,X) are solutions of (|23J) that decay in 5*0, we see 
that these two are linearly dependent. So one is a constant multiple of the other. Moreover, 
from we see that b}~(a, A) = bk(a, A) for all k G N where we used bk(a, A) instead of bk 
to indicate its dependence on a and A. Also we have F(z, a, A) = F(z,a, A) in Proposition 
f|. Hence the entire functions g(z) and f{z,a, A) along with their first derivatives satisfy the 
same asymptotic expressions in Proposition § (ii), so we conclude that 



(24) g(z) = f(z,a,X) 

by Proposition |] (iii). Next substituting z = in (p4|) gives the first equation in (0). Also 
we differentiate (p4|) with respect to z and substitute z = to get the second equation in 



(0). For just note that G(a) = G _1 (a). □ 

Next we want infinite product representations of /(0, a, A) and /'(0, a, A), with respect to A. 
But first, we recall the definition of order of an entire function, which will be needed in the 
proof of the next lemma. Let M(r,g) = max{|(yf(re je )| : < 9 < 2tt} for r > 0. Then the 
order of an entire function g is 

log log M(r,g) 
lim sup . 

r-,00 log r 
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If for some positive real numbers a, c±, C2, we have M(r, g) < c\ exp^r ] for all large r, then 
the order of g is finite and less than or equal to a. 

Lemma 8. Suppose m > 3. The functions X ^ /(0, a, A) and A 1— > /'(0, a, A) nave infinitely 
many zeros Ej and Ej, respectively. They admit the following infinite product representations 
for each fixed a = (ai, a 2 , . . . , a m -i) G C m_1 : 

00 / A A 

/(0, a, A) = D X n ° J^J I 1 j for some D 6 C and nonnegative integer n , 

3=1 ^ j ' 
00 / A A 

/'(0, a, A) = _DiA ni 11(1 7 I /or some Z)i 6 C and nonnegative integer n\. 

Moreover, these infinite products converge absolutely. 

Proof. If we show that both /(0, a, A) and /'(0, a, A) have orders (with respect to A) strictly 
less than one, then this lemma is a consequence of the Hadamard factorization theorem 
(see, for example, Theorem 14.2.6 on page 199 of [15[]). So we will show that /(0, a, A) and 
/'(0, a, A) have orders strictly less than one. 

From the equations (|10[) and (PI), we see that except for 7r — <5 < arg A < it + 5, small 5 > 0, 
both |/(0, a, A) I and |/'(0,a, A)| are bounded by exppi^jAI 1 ^?] for large |A|. So to show that 
they have the orders strictly less than one, it suffices to show that for tt — 5 < arg A < it + 5, 
small 5 > 0, they are bounded by d\ expf^lA] 1 ^] for some d\ > and c?2 > 0. 



From (|i~4] ) with z = 0, one can see that 

|/(0, G-^a), A)| < \C(a,u 2 X)f(0,a,uj 2 X)\ + C{a,uj 2 X)f{0,G{a),uj- 4 X) . 

In this inequality, let A lie in the region | argA — 7r| < 5. Then since a; 2 A and uj~ 4 X are not 
in I argA — 7r| < 5, we can use the asymptotic expression ( P~D[ ) to get that for all large |A|, 

|/(0, ^^(a), A)| < [|C(a,w 2 A)| + C(a,u 2 X) j exppKlAI^ 2 ]. 

We know that C(a,co ,2 A) depends only on a by Lemma [5| Also the equations (29.4) and 
(29.7) imply that for fixed a, 

\C(a,u 2 X)\ < d3exp[2if|A|^], for some d 3 > 0. 

Thus we see that for each a, 

1/(0, G-\a),X)\ < diexp[4|A|^r], for some d\ > 0. 



m+2 
2m ' 



Hence the order of /(0, a, A) with respect to the A-variable is less than or equal to 
Hence by combining this with (|I0"D, we conclude that the order of /(0, a, A) is 2 y^, which is 
strictly less than one. 

Next we differentiate ( |HD with respect to the ^-variable and set z = 0, then similarly 
using (|n|), we can conclude that the order of /'(0, a, A) is □ 
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3. Proof of Theorem 
When m — 2, the equation ([!]) is a translation of the harmonic oscillator. So there is 



nothing new here. We mention that since z 2 + a±iz 



■1 > 



the eigenvalues for the 



potential z 2 + a\iz are Ik + 1 + ^ > 0. 

Suppose m > 3 and suppose that A G C is an eigenvalue of the eigenproblem ([[]), then by 
Lemma ^| we have C(a, A) = 0. Then from ( jl9|) and (^) along with (|22|), we have 

= /(o, G _1 (a), w m A) - e^°/(0, G _1 (a), cu m A), 

= w/'(0, G _1 (a), cj m A) - e^u^f'iO, (^(a), u; m A). 

Since the non-constant function /(£, G _1 (a), c<j m A) solves a linear second order ordinary 
differential equation, both /(0, G _1 (a), co> m A) and f'(0,G~ 1 (a),u m X) cannot be zero at the 
same time; otherwise, f(z, G _1 (a), u) m \) = 0. 

Suppose that /(0, G^ 1 (a),uj m X) ^ 0. Then from Lemma |S] we have 



o^^o^Arn 

3=1 



1 - 



E; 



i=i 



u; m A 



Then by equating the absolute values of the two sides of the equation (and using u m+2 = 1), 
we have 

u 2 Ei - A 



(25) 



n 

i=i 



1. 



uj 2 E j - A 

Likewise, when f'(0, G~ 1 (a), uo m X) ^ 0, we get the following. 



(26) 



n 



u 2 E' - A 



^E] - A 



1. 



We mention that uj 2 Ej and ^E 1 -, lie in the open lower half-plane for some j, j'. From 
Lemma [8] we know that f(0, G _1 (a), E) and f'(0, G~ l (a), E) have infinitely many zeros E*. 
And (|10D and (0) imply that the zeros £"* near infinity lie near the negative real axis. Thus 
certainly \muj 2 Ej < and \muj 2 E'j, < for some j, j'. 

Below we will show that the hypotheses on the signs of the coefficients a 1; a 2 , . . . , a m _i of 
P force all the uj 2 Ej and u 2 E'^ to lie in the closed lower half-plane, which implies either 

Vj G N, if ImA > 0, or 
Vj G N, if ImA < 0, 

since A and A are reflections of each other with respect to the real axis. If (|25|) holds, 
then (|3) implies \oo 2 Ej - A| = \uo 2 Ej - A~| for all j G N. If © holds, then fl27|) implies 
|o; 2 ^ — A I = l^E'j — A | for all j G N. Since Im u 2 Ej < for some j and Imu 2 ^ < 
for some j, and since A and A are reflections of each other with respect to the real axis, we 
deduce in either case that A = A and hence A is real. 



\uo 2 Ej - A 


> 


\oo 2 Ej - A 


and 




- A| 


> 


\u 2 E> 


- A| 


\uo 2 Ej - A 


< 


\u 2 Ej - A 


and 


\. 2 E< 


- A| 


< 


\u 2 E> 


- A| 
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So our next task is to show that all the w 2 Ej and oj 2 E'- lie in the closed lower half-plane. 
Suppose that for some E* G C, 



(28) 



either /(0, G _1 (a), Ej = or /'(0, G _1 (a), E») = 0. 



That is, either .E* = Ej or E* = E'- for some j G N. We know that v(z) = f(z,G 1 (a), E*) 
solves 

(29) 



-v"(z) + 



m—l 



+ a k u k z m - k 



k=l 



v(z) = —E*v(z), 



where a k G JR., k — 1, 2, . . . , to — 1, and by (|28|), i> satisfies either Dirichlet (i?* = Ej) or 
Neumann boundary = Ej) condition at 0, and Dirichlet condition at oo + Oi. We call 
(p9|) with these boundary conditions the "associated" eigenvalue problem. We aim to show 
all the eigenvalues E* have Im (uj 2 E*) < 0. 

Let g(r) = v{re %e ) with 6 fixed, \6\ < We then multiply ( |2"9"D by u) 2 e~ 2ie g{r) and 

integrate over < r < oo to get 

m— 1 



2 -2i0 



(30) -u/e 



-CJ 2 K 



g"{r)g(r) dr + 
|<?(r)| 2 dr. 



u 2 e mie r m 



+ {a k u k+2 e^- k ^ e r m - k ) 



k=l 



\g(r)\ 2 dr 



Since f(z, G _1 (a), E 1 *) decays to zero exponentially in Sq, we know the integrability of every 
term in ( |30"D for each \8\ < ^jy. Next we integrate the first term by parts, using g(0) = 
or <7'(0) = by (^), so the boundary term vanishes. And then taking the imaginary part of 
the resulting equation gives 



sin 

m— 1 



47T 



TO + 2 



20^ jf \g'(r) | 2 dr + sin f m0 + 



47T 



TO + 2 / ,/ 



r m |#(r)| 2 dr 



+ £ 



afc sm 



fc=i 



2(A; + 2)tt 
to + 2 



(31) 



+ (m - k)9^j jf 

/■oo 

-Im (w 2 ^) / |#(r)| 2 dr. 



r m - fc |#(r)| 2 dr 



Recall our hypothesis that (j — k)a k > for all 1 < k < m — 1 for some 1 < j < ~. We 
want to prove the reality of the eigenvalues by showing that Imu 2 £» < for all the E*. To 
this end, we will divide the proof into two cases; Case I, when 1 < j < — and m > 5, or 
when j = 1 and m = 3, 4; and Case //, when j = 2 and to = 4. 

Case /: when 1 < j < y and m > 5, or when j = 1 and to = 3, 4. We choose 6 in ([H]) by 

(to - 2j - 2)vr 



(32) 







(to — j)(m + 2) 
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where the motivation for this choice will be fairly clear later in the proof. Notice here that 
\6\ < ^2 as required. Then 

4vr „ 27r , , / 47r \ 

< 26 = < 7T, and hence sin 26 > 0, and 

m + 2 m — j \m + 2 J 

An (m - 2j)n , , / n An \ 

< m6 H = < n, and hence sin mv H > 0. 

m + 2 m — j \ m + 2 J 

(Clearly these inequalities use that j < — .) Also we see that 

2(fc + 2)7r + (m-fc)fl = ( 1 _tll\ vr, forallfc, 
m + 2 \ m — j J 

and so 

sin (l — — J 7r has the same sign as (j — k) for all 1 < k < m — 1. 

\ m-jj 

Among other things, this is why we choose the 6 as above. So from ( [31]) and the hypothesis 
(j — fc)afc > 0, we conclude that Im (u 2 E*) < 0. This proves that the eigenvalue A is real. 

Case II: when j = 2 and m = A. The reason we separate this case from Case I is that in 
this case, \6\ = ^ = ^r^, whereas our argument needed \6\ < in order to get integrability 
for terms in (^3). So we modify the proof as follows. For e > small, we multiply (|30|) by 
e~ 2le and set 6 = — | + e. Then integrating the first term by parts and taking the imaginary 
part of the resulting equation give 



roo poo 

sin(4e) / |g'(r)| 2 dr + sin(2e) / r 4 \g(r)\ 2 dr 
Jo Jo 



(33) + £a fc ffln^ + (2-A;)eJ jf r 4 "^(r)| 2 dr 

POO 

= -Im (u 2 E,e- 2ie ) / \g(r)\ 2 dr. 
Jo 

Clearly sin (If + (2 - k)t) has the same sign as (2 — k) for 1 < k < 3. Using the hypothesis 
(2 — k)a,k > for all 1 < k < 3, we have that the left side of ([33]) is nonnegative and so 

Imu 2 B t e" 2,e < 0. 

Thus by sending e to zero, we get 

Imw 2 £ t < 0, 

which proves the reality for the case of j = 2 and m = A. 

Therefore, the eigenvalues of (P are all real under the hypotheses on the a^'s given in the 
statement of this theorem. 

We must still prove the positivity of the eigenvalues. Suppose u is an eigenfunction of 
(H|) with an eigenvalue A 6 1, and suppose ak's satisfy the hypotheses of the theorem. Let 
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v(z) = u{—iz). Then we have the equation ([8]) with the boundary conditions that v decays 
in 



S_! uS 1 = {ze 



7T 



m + 2 



< are; z\ < 



3tt 



m + 2 



}• 



Since A and all a^'s are real, one can see that v(z) satisfies the same equation and decays in 
S-i U Si. Then since the eigenvalues are simple, v{z) and v(z) must be linearly dependent, 
and hence v(z) = cv(z) for some c £ C. Since |t> (z)| and agree on the real line, we see 

that |c| = 1 and so \v(z)\ = \v{z)\ for all z £ C. That is, |f (x + is even in y. From this 
we have that 

(34) 



d 

= a~\ v ( x + ^)l 2 

oy 

J.8 



y=0 



-21m (v'{x)v (x) j , for all x £ R. 



Next we let /i(r) = v(re ). By substituting into the differential equation @, then multi- 



plying by h(r) and integrating, we get 



h"(r)h(r))dr + 



m—l 



ake ( m -k+2)i8 r m-k 



k=l 



\h{r)\ 2 dr 



poo 

-Xe 2ie \h(r)\ 2 dr, for <9 

Jo m -\- 2 



< 



3tt 



m + 2 

Integrating the first term by parts and using h'(0) = e j6 V(0), one can get 

m—l 

|//(r)| 2 dr + 



i/(0)u(0) + e 



(m+l)i6^,m _|_ ^ ^ ^^{m—k+VjiO^m—k 



k=l 



\h{r)\ 2 dr 



POO 

-Xe ie / \h(r)\ 2 dr, for 



7T 



m + 2 



< < 



3tt 



m + 2 



Taking the imaginary part and using (|34|) at s = 0, we have 



sin 



m—l 



r sin m 



+ 1)0 + a kr m ~ k sin(m - k + 1)0 



k=l 



\h\ 2 dr 



Asin0 / \h\ 2 dr, for all 



7T 



(35) , , , 

m + 2 

(Here again we used that A is real.) We choose 





< 6 < 



3tt 



m + 2 



7T 



so that 



as required, and 



7T 



< 



7T 



m — j + 1 ' 
2?r 



m + 2 m + 1 



< < 



3tt 

< < 7T 



sin(m — k + 1)0 = sin 



m + 1 m + 2 
m — k + 1 



-7T 



m — j + 1 

has the same sign as (k — j), for all /c. Since (k — j)a,k < 0, sin0 > and sin(m+ 1)0 < 0, we 
see that the left-hand side of (|35|) is positive, and hence so is the right-hand side. Therefore, 
the real number A must be positive. This completes the proof of Theorem 0. 
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Remarks. 

1. The idea of using the infinite product in (|25|) to prove reality of the eigenvalues is due 



to Dorey et al.||ll]|. But their potentials are much simpler and the E* are all negative real 
in their situation, so that fl2~7|) is immediate. Here is not. 

2. The ideas above for proving positivity of the eigenvalues are similar to those used earlier 



by the author in |TJ| . 

3. We note that the hypotheses assumed in Theorem ^| on the coefficients of P are sufficient 
for real eigenvalues, but not necessary, for at least two reasons. Let Q(z) = — \z m + P(z)]. 
Then first, the problem (fl|) with the potential Q(iz) = — [(iz) 3 — (iz) 2 ] is covered by Theorem 
H] while the problem with Q±(iz) = —[(iz) 3 + 2(iz) 2 + iz] is not. However, Q(z + 1) = Qi(z) 
and so the potential Qi(iz) produces positive real eigenvalues only. For general cases, for a 
real polynomial Q if the problem ([!]) with the potential Q(iz) has positive real eigenvalues 
A only, then the problem with the potential Q(iz + c) — Q(c) for some real c G IR has 
eigenvalues A — Q(c) which are all real. Second, in the proof of Theorem || in order to ensure 
that Imcu 2 ^ < 0, we insisted that each and every term on the left-hand side of (|31]) has a 
single sign, and it is clear from Section f| below that this is not necessary. 

4. Extensions of Theorem || 

In this section, we study two particular classes of polynomial potentials to illustrate dif- 
ferent methods for sharpening Theorem |2|. 

Theorem 9. Let m > 4 and suppose a < 0, 7 < 0. Suppose that an entire function u along 
with A G C solves the equation (fjj with P(z) = az 3 + j3z 2 + jz. Then the eigenvalue A is 
positive real, provided that 



(36) /3 < v^7y3 - tan 2 (J^j . 

The eigenvalue A is also positive real provided that A el and 



tan 2 

(37) (3 < 4v / 2 v ^7- 



3-tan 2 (^-) 

V m+l / 

Remarks. 

1. Note that for a, [3, 7 < 0, we have A > by Theorem |2|. The point of Theorem^, then, 
is that if a, 7 < then we can allow some values of (3 > 0. 

2. The right-hand side of fl3"6] ) is less than that of ( |3"TD as we show at the end of the proof. 

Proof of Theorem |^. Since the theorem for (3 < is contained in Theorem [2|, it suffices to 
show the claims of the theorem hold under the hypotheses (|36|) and (|37|) with (3 replaced by 
|/3 1 . In proving this we will closely follow the proof of Theorem |2|. 
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As in the proof of Theorem |2|, in order to prove the reality of the eigenvalues we show that 
Im (o; 2 ^) < for all E* G C satisfying fl28|). In this case, we see that (|3T|) becomes 



oo 

sin (20) / \g'(r)\ 2 dr — sin(m(p) / r m \g(r)\ 2 dr 
Jo Jo 

/■OO 

(38) — / [ar 2 sin (30) + (3r sin (20) + 7 sin 0] r|g(r)| 2 dr 

/■oo 

= -Im (cj 2 £*) / \g(r)\ 2 dr, 
Jo 

where = — 9. (See the proof of Theorem § for the definition of g{r).) 

Recall that the positivity of the left-hand side of (^) implies A Gl. Since \9\ < ^r^i we 
get < < Then since we are trying to show the left-hand side is positive under 

certain conditions on the coefficients, we restrict to — < < — if m > 5, and ? < < £ 
if m = 4, so that sin(m0) < in the second term above. (Note that when m = 4, = ~, 
we have = — | for which some terms in (|38|) are not integrable.) We further want the 
discriminant of the quadratic [ar 2 sin (30) + j3r sin (20) + 7 sin 0] to satisfy 

I3 2 sin 2 (20) -4ct7sin(30) sin0 < 0, 

so that the quadratic expression has a single sign. That is, we want 

(oc\\ (V ^ a sin(30)sin0 , 2 . 

(39) /3 < 4«7 — . 2 = 07 (3 - tan 0) . 

sm (20) 



So in order to have a positive right-hand side in (|3~9"D, since a < and 7 < 0, we need 
^ e [m' f) ^ tm' ~] ^ or wn i c h (3 — tan 2 0) is positive. Since (3 — tan 2 0) is decreasing, 
to maximize the right-hand side of (|39"D, we choose = — . Hence as we remarked at the 
beginning of the proof, this proves the reality of the eigenvalue under (j36|). 

Similarly, in order to prove the positivity of the eigenvalues, suppose A G E and use (|35|). 
Let ft.(r) = v(re ld ) = u(—ire 10 ). Then (pop becomes 



/■OO /■oo 

sin0 / |/i'| 2 dr-sin(m + l)0 / r m \h\ 2 dr 
Jo Jo 

/■oo 

(40) - / [ar 2 sin (49) + (3r sin (39) +7 sin (29)] r\g(r)\ 2 dr 

Jo 



7T 37T 

A sin / \h\ 2 dr, for all < 9 < ' " 







m + 2 m + 2 



Then we restrict G [— tt, - 2 tt1 so that sin(m + 1)0 < in the second term above. We also 

Lm+l' m+lJ \ ' — 

want the discriminant of the quadratic [ar 2 sin (49) + (3r sin (30) + 7 sin (20)] is nonnegative, 
so that the quadratic expression has a single sign. That is, 



/? 2 sin 2 (30) -4a7sin(40)sin(20)sin 2 (30) < 0. 
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Hence we have 



sin(4fl) sin {29) 1 - tan 2 9 



(41) p A < 4a7 — v y 9/ „; 7 = 32a 7 - 

V ; sin 2 (3fl) (3-tan 2 #) 

One can check kzten?e ■ decreasing on [— tt, f) . Also we want to have (1 — tan 2 9) > 

(3-tan 2 6») to L m+l ' 3/ V / — 

so that the right-hand side of (|4*T| ) is nonnegative (that is, T^pj < 9 < |.) Then it is not 
difficult to see that 9 = maximizes the right-hand side of fl4"Ir). Also with 9 = — tt, the 

m+l ° \L_J/ m+l ' 

left-hand side of ((4"0|) is positive and hence A > 0. So with help of Theorem |2|, we conclude 
that all real eigenvalues are positive under the hypothesis (|37|) . 

Still we must show that eigenvalues are positive under (|3~6"D . We will do this by showing 
that the hypothesis ([36]) implies fl37l) . That is, we will show 



o/7T\ l-tan 2 (^ T ) 

(42) 3 - tan 2 ( - ) < 32 Wi; 9 , for all m > 4. 

U; (S-tan 2 ^)) 2 ' 



Since ^ is decreasing and positive for 9 G [0, |) , the right-hand side of (^) is an 

increasing function of m > 5 and hence greater than or equal to the value at m = 5 which 



is 3. So ( [42]) holds for m > 5 since its left-hand side is less than 3. And for m = 4, one just 
check ([42"D directly. This completes the proof. □ 



Remark. Above we have chosen P(z) = az 3 + (3z 2 + 72; for simplicity. One should note 
that the above argument works for real polynomials of the type P(z) = az n+k + (3z n + / yz n ~ k 
for some positive integers n > k. 

The previous theorem handled a < 0. Similarly, we get the following for a > when 
m = 4, 5, 6. 

Theorem 10. Let m = 4, 5, or 6 and let a > 0, 7 < 0. Suppose A is an eigenvalue of ([![) 
P(^) = az 3 + /?2 2 + 72 /or some /3 G R. Tnen tne eigenvalue is positive real, provided 

that 



(43) $ < Wtan 2 — - 3 



/?< v^Rjtan 2 0£) 



Tne eigenvalue A zs a/so positive real provided that A G R and 

00 if m = 4, 5, 

(44) /3 < 



'l-tan 



(¥) 



Remark. In this theorem we restrict m to m = 4, 5, 6 for reasons explained in the proof 
below. 



In ([43D, by convention we take a/|«7| y tan 2 (^) — 3 = +00 when m = 4, so that ( [43"D 
just says /3 G R, in that case, as in Theorem |2[. Note that by (|44]), all real eigenvalues are 
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positive when m = 5, a > 0, (3 £ R, 7 < 0. But there could perhaps be some non-real 
eigenvalues. 



We mention that the case (3 < or m = 4 of Theorem |10| is contained in Theorem |2], as 

are new. 



it is explained in the proof of Theorem 9L The other cases of Theorem 10 



Proof of Theorem |7^. Since the case (3 < is known already, it is enough to prove the 
theorem under the hypotheses (|43"D and ([44] ) with (3 replaced by \(3\. 

The proof below will be very much similar to that of the previous theorem. So we will 
refer equations to those in the proof of the previous theorem. 

Since the case m = 6 in ( fi3|) says (3 < 0, this case is contained in Theorem [| So for ( f4~3D 
we can assume m = 5. Again we use ( pi]) with P(-z) = az 3 + /3z 2 + 72;. Then as we did in 
the proof of Theorem |9|, we can get fl3~8|) , where we want - < < — so that sin(m0) < in 
the second term. But this time since a > 0, we want sin(30) < and want the discriminant 
of the quadratic [ar 2 sin (30) + (3r sin (20) + 7 sin 0] to be non-positive. Then again we have 
(|55D. Obviously we want a nonnegative right-hand side in ([39]), and since a > and 7 < 0, 



we need (3 — tan 2 0) < 0. This means | < < | as well as ^ < < — . 

Then in order to maximize the right-hand side of (p9|) we choose = ^ for m = 5. This 
proves the reality of the eigenvalue since (3 < is covered by Theorem |2|. 



Similarly, in order to prove the positivity of the eigenvalues under the hypothesis (|44j) , 
suppose m = 5 or 6, let A e M and use fl35D . Then like before we get (^l|) where we want 
Q e [^1, ^fr] so that sin(m + 1)^ < in the second term in (0). We further want 
the right-hand side of (f|l]) to be nonnegative, and hence want > |. For m = 5 since 
f < f < ^j, we get /3 < 00 from (0). And for m = 6 since f < ^ = ^ < f , and 
since , 1 ~ tan e A < is decreasing for 9 £ |) , we choose # = ^ in order to maximize 

(3— tan^ 0) L 4 o / r 

the right-hand side of (|4T|). Thus along with Theorem ^| for /3 < 0, we conclude all real 



eigenvalues are positive under the hypotheses AeR and (|44|). 

Finally, it is not difficult to see that fl43|) implies (|44|) , and hence we get the positivity of 
the eigenvalue under ([43]) as before. This completes the proof. □ 



The second method for sharpening Theorem [^ is to make use of the J \g'(r)\ 2 dr term, by 
means of the harmonic oscillator inequality. Below we will prove that eigenvalues A are 
positive real if a < y. Note that Dorey et al. already prove this, and they also show 
that eigenvalues are real if a < y + 2. 

Theorem 11. Let m > 4 be an even integer. Suppose that an entire function u along with 
A G C solves 

d 2 



(45) 



dz 2 



- (iz) m - aiiz)^ 1 



u(z) = Xu(z) 
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with the boundary conditions that u decays to zero as z tends to infinity along the rays 



&rgz 



■5 ± 



2 IT 



,, ^_ (i , 2 . Then the eigenvalue A is positive real if a < y. If a = ^, then all 
eigenvalues are positive real except the smallest one, which is zero and has the corresponding 
eigenfunction 



u Q (z 



Note that the function v(z) = exp 



cxp 



m+2 



m + 2 



m + 2 



Z 2 



solves — v"(z) + z m v(z) 



-Z 2 



"^(z) and 



decays along both ends of the real axis. This type of problems was studied by Bender and 
Wang §. 

Proof. The outline of the proof is similar to those of the proofs of Theorems ^ and [10 above. 
But this time we will make use of the f \g'\ 2 dr term via the harmonic oscillator inequality. 
Like before, we examine the equation (|3l|) , for the choice P{iz) = aiiz)^ -1 . Then we have 
An __\ f°° . / „ An 



sin 



— a sin 



m + 2 

An 



m + 2 



20 I I \g'(r)\ 2 dr + sin (m9 + 

J Jo ' 1 )l V ™ + 2 

r ^ 1 1 c?(t) | 2 dr 



r m \g(r)\ 2 dr 



o 

oo 



(46) 



POO 

-Im (aAE*) / \g(r)\ 2 dr 
Jo 



where we refer to the proof of Theorem |2| for the definition of g(r). Now we use the harmonic 
oscillator inequality on the first two terms above so that we have that for \9\ < -^p^ with 

s in H + ^) > °> 

An \ f 00 „ . / „ An 



sin 



m + 2 



2d ) So \ 9 '( r ^ 2dr + sin ( m0 + 



> 



> 



sm 



An 



m 



29 sin m9 + 



An 



m + 2 



sm 



47T 



m 



29 sin m0 + 



47T 



m + 2 



ml [An 

— \ / sin 

2 V 1 m + 2 



29 sin m# 



47T 



m + 2/ j 

by parts. We then combine this with ( pf) to get that if 



m + 2 
2r~\g'(r)g(r) \ dr 

f°° m d 

I r ~ \n( r \\ 2 a] r 
Jo dr^ y 

r^~ l \g(r) \ 2 dr, 



r m \g(r)\ 2 dr 



■ m I . ( An 

(47) — i sin 

2V \m + 2 



29 sin m6» + 



An 



> a sin 



47r m - 2 „ 

+ —^0 



m + 2 J \m + 2 

then Im (uj 2 E*) < 0, which then proves the reality of the eigenvalues like in the proof of 
Theorem |2]. Next we examine the condition (f47|) and find with a little effort that 9 = is 
the best choice to get the best bound for a out of (|47D . That is, a < y. 
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Similarly, in order to prove the positivity and non- negativity of the eigenvalues, we use 



). Let h(r) = u(—ire ). Since A G M, one can get the following from (| 

.9 f \h'\ 2 dr-sm(m+l)9 [ r m \h(r)\ 2 dr - asm C^B) [ 
Jo Jo ^ 2 / Jo 



sm i 



f 2 



-1 



h\ 2 dr 



A sin# 



\h\ 2 dr, 



provided 



7T 



< e < 



3n 



Then since sin(m + 1)9 < for 
to the first two terms above to get 



m+1 



< 9 < 



m + 2 m + 2 

we apply the harmonic oscillator inequality 



> 



> 



A sin# 



rn 
12 
m 
~2 



\h\ 2 dr 



m 



sin (9 |sin(m + 1)9\ — asm -^9 



™ 1 \h\ 2 dr 



where we have chosen 9 



a 



2~ 



T 2 



'^hfdr, 



. (One can check with a little effort that 9 



2- 



is the best 



choice for this argument.) Since sin > 0, we see A > when a < y and A > when a = y . 

When a = y, we see that for m even, Wo(— iz) solves —v"(z) + (z m + y^" 1 ) = 
with properties that uo(-iz) decays in S-i U Si and blows up in Sq. So it satisfies the proper 
boundary conditions of (§) and hence u (z) is, in fact, the eigenfunction of (^5|). Hence, all 
eigenvalues are positive except the smallest eigenvalue zero since eigenvalues are simple by 
Proposition |I|. This completes the proof. □ 



Remark. Note that arguments similar to Theorem [IT] work, for example, for some poly- 
nomial potentials with P(z) = (3z 2 + 72 when m > 7. Note also that to prove the reality 
and positivity of the eigenvalues, it is enough to show the left-hand sides of (31) and (S5j) 



are positive, respectively. For specific potentials, this might be achieved by other kinds of 
estimates. 

Remark. The methods of proving the theorems in this section show how to sharpen The- 
orem ^| to problems with potentials "almost the same" as those of Theorem |^. But for 
Theorem [11], the proof of the reality can be shortened as follows. 

Suppose that an analytic function v(z) = f(z, G _1 (a), E*) along with 6 C solves 



-v"(z) + (z m - azf-^viz) 



-E*v(z), with v'(0)v(0) = 0, v(+oc + Oi) = 0. 



This eigenproblem on the positive real axis is self-adjoint, and so E m G M. Precisely, mul- 
tiplying both sides by v(z), integrating over the positive real axis, and integrating the first 
term of the resulting equation by parts give 



\v'(x)\ 2 



dx + I x m \v(x)\ 2 dx — a / x 2 1 \v(x)\ 2 dx 
o Jo 



-E ± 



\v (x) \ 2 dx. 
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(Hence E* el.) Then one uses the harmonic oscillator inequality on the first two terms to 
have 

(lT ~ °) J x ^~ 1 \ v i x )\ 2 dx < ~ E * J \v{x)\ 2 dx. 

So if y > a then E* < which implies Im {uj 2 E*) < 0. Thus the eigenvalue A is real, as 
before. 

5. Conclusions 

In this paper we have proved that a family of one dimensional Schrodinger equations with 
PT-symmetric polynomial potential — [(iz) m + a^iz)"" 1 ^ 1 + a2(iz) m ~ 2 + ■ • • + a m _i(iz)] has 
all positive real eigenvalues, provided that (j — k)ak > for all k, for some 1 < j 1 < y. In 
particular, this result implies the original Bessis and Zinn- Justin conjecture for the potential 

iaz 3 + z 2 . 

One would like to further extend the proof of the reality and positivity to a larger class 
of PT-symmetric potentials. For example, can we get a similar conclusion for j > y? 
Also an interesting question will be how much the reality and positivity of the eigenvalues 
depend on the boundary conditions: one certainly has some restrictions on choosing the 
boundary conditions (see, for example jp). Our boundary conditions allow only one blowing 
up Stokes sector between the two decaying sectors near the negative imaginary axis. It will 
be also interesting to consider three or more (an odd number of) blowing up sectors between 
the two decaying sectors on which we impose the boundary conditions, with the decaying 
sectors being symmetric with respect to the imaginary axis. Also the problem with the 
potentials +{iz) m — P{iz) (whose leading term has the opposite sign to those in the class of 
the problems studied in this paper) would be interesting too, in which case we impose the 
boundary conditions to allow an even number of blowing up sectors between the decaying 
sectors. Also it should be possible to apply the methods of this paper to some rational 
potentials, too. 

One big question needing to be answered is to determine the span of the set of the eigen- 
functions. For Sturm-Liouville problems, we know that zeros of eigenfunctions interlace, 
which seems to play an important role in the completeness of the set of the eigenfunctions. 
Numerical work of Bender et al. || shows some intriguing interlacing properties of the ze- 
ros of the eigenfunctions for some PT-symmetric oscillators, too. So understanding these 
interlacing properties of the zeros might lead us to progress. But yet, it seems there are a 
lot more questions than answers in this direction. 

Also, one would like to study similar problems in higher dimensions. 
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